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Abstract 



We have investigated the far-infrared dipole modes of a quantum antidot 
submitted to a perpendicularly applied magnetic field B. The ground state of 
the antidot is described within local spin-density functional theory, and the 
spectrum within time-dependent local spin-density functional theory. The 
results are compared with those corresponding to a quantum dot of similar 
electronic surface density. The method is able to reproduce two of the more 
salient experimental features, namely that main bulk and edge modes have the 
same circular polarization, and that the negative B dispersion edge branch 
oscillates, having minima at the B values corresponding to fully occupied 
Landau levels. It fails, however, to yield the unique feature of short-period 
antidot lattices that the energy of the edge magnetoplasmon approaches the 
cyclotron frequency for small B. The existence of anticyclotron polarized 
bulk modes is discussed, and a detailed account of the dipole spin mode is 
presented. 

PACS 73.20.Dx, 73.20.Mf, 78.20.Bh 
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I. INTRODUCTION 



The far-infrared (FIR) spectroscopy of laterally confined superlattices in the two di- 
mensional electron gas (2DEG) made of holes surrounded by electrons, called antidots, has 
uncoveredEhi a very peculiar behavior of the collective spectrum of these systems. It has 
been foundllli that it consists of a high frequency branch which starts either with a nega- 
tive magnetic field (B) dispersion in the case of high electronic surface density n^, or with 
a rather flat B dispersion at low Ug, that eventually converges to the cyclotron energy at 
high 5's, plus a low frequency branch which at high 5's corresponds to the usual edge 
magnetoplasmon, but that approaches the cyclotron frequency for small magnetic fields. 
Anticrossing of the modes appears as B increases. It is worth to recall that the mentioned 
behavior of the high frequency branch is also a feature of the FIR spectrum of quantum 
rings§li. 

A study employing circularly polarized radiation! has shown that both edge and bulk 
magnetoplasmons in antidots exhibit the same circular polarization, in contradistinction 
with what is found in quantum dotsB'i. Interestingly, a very recent experiment! has detected 
weak bulk modes anticyclotron-like polarized, whose existence had been predicted some time 
agoSini. 

Another difference with dots appeared after a careful analysis of the edge 
magnetoplasmon0: the frequency of this mode shows a conspicuous oscillation with B 
that has maxima at fully occupied Landau levels (even filling factors u) in the case of dots, 
and at half-filled Landau levels (odd i^'s) in the case of antidots. 

Several theoretical descriptions of FIR modes in antidots have been given in the past. 
Some are based on hydrodynamicaSEl or classical electrodynamic^H models, which either 
do not take into account the periodicity of the antidot arravu, or which take it into account 
either in a circularly symmetric Wigner-Seitz approximationt°l, or incorporate the parameters 
of the experimental short-period antidot lattices without further approximationsl^. 

Despite that all these models have somehow succesfuUy reproduced the FIR absorption 
in antidots, the interest in achieving a more microscopic description of the ground state 
(gs) and FIR response still remains. Quantum mechanical (Hartree) models have been 
set up which take into account the periodicity of the lattice using confining potentials of 
different complexitylii^. The difficulties inherent to the handling of many single particle 
(sp) wave functions and realistic confining potentials which refiect the actual antidot lattice 
have hampered these microscopic methods to achieve a quantitative description of the FIR 
absorption in antidots, although some features of the process are qualitatively describecllil'0. 

Recently, we have applied a local spin-density functional (LSDFT) method to calculate 
the structure of isolate antidots at B = 0, together with a sum rule approach to describe 
their FIR spectrum0. Our aim here is to extend these calculations to the case of finite 
magnetic fields and to present an account of the gs and dipole response of an antidot within 
the frame of LSDFT and its time-dependent generalization (TDLSDFT). We recall that 
these formalisms incorporate the spin degree of freedom and allow to selfconsistently include 
exchange and correlation effects in the study of the gs and the response of the system, whose 
importance in the description of 2D electronic structures is nowadays well established. 

Rather than isolated antidots, the systems described in this work should be considered 
as representing long-period antidot lattices. Although no much differences are expected to 
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appear between the ground state structure of an isolated antidot and that of an antidot 
in a long-period array, apart from the obvious changes at the border of the unit cell, bulk 
magnetoplasmons are qualitatively different in both case^. 



II. THE LSDFT AND TDLSDFT APPROACHES 

We have modeled an antidot of radius R in a. 2DEG of surface density Ug by a positive 
jellium background of density nj{r) = ns0(r — R). The gs of the antidot is obtained 
solving the Kohn-Sham (KS) equations. The problem is simplified by the imposed circular 
symmetry, which allows one to write the single particle (sp) wave functions as (pniaif-iO) = 
Uni(7{f)€-~^^^ with / = 0, ±1, ±2, being —I the sp orbital angular momentum. 

We have used effective atomic units defined hy h = e^/e = m =1, where e is the dielectric 
constant, and m the electron effective mass. In units of the bare electron mass irie one has 
m = m*me. In this system of units, the length unit is the effective Bohr radius Og = aoe/m*, 
and the energy unit is the effective Hartree H* = Hm*/e^. In the numerical applications 
we have always considered GaAs, for which we have taken e = 12.4, m* = 0.067, and 
g* = -0.44. This yields ~ 97.94 A and H* ~ 11.86 meV ~ 95.6 cm'^ The Bohr 
magneton is fiB = fie/2meC, and the cyclotron frequency is Uc = eB/mc. 

The radial KS equations read 

1 f (f Id /^\ tUc, 1 2 2 / N 

-2[d^^^rd-r--^)-Y^ + r^' 

(1) 



where rjcr = +1(— 1) for cr =1 (|), = J d'r'n{f')/\f—f'\ is the Hartree potential, and 
yxc _ Qg^^iji^fji^ I Qfi\^^^ and W^^ = dSxc{n,m) /dm\gs are the variations of the exchange- 
correlation energy density Sxc{n,m) written in terms of the electron density n(r) and of 
the local spin magnetization m(r) = n^(r) — n^{r) taken at the ground state. £xc{n,m) 
has been built from the 2DEG calculations of Tanatar and CeperleyEl at zero and full spin 
polarization, using the Von Barth and HedinS interpolation formula adapted to 2D. The 
jellium potential V~^{r) is analytical: 

+ r i?ooE(r/i?o,) - i?E(r/i?) r<R 

V [T) ^^-^ j/^^E(r/i?oo)-rE(i?/r) + r[l-(i?/r)2]K(i?/r) r>/2, ^ > 

where K and E are the complete elliptic integrals of first and second kind, respectively^, 
and Roo represents a large r value. In practice, it is the largest-r used in the gs calculation. 
The KS differential equations have been solved without expanding the sp wave functions 
in a necessarily truncated basis of Landau orbitals. Our iterative method works for weak 
and strong B fields as well, for which the effective potential is very different. It has the 
advantage of avoiding to study how the results depend on the size of the basis. As in Ref. 
we have worked at a small but finite temperature (0.1 K). 

Physically acceptable solutions have to be regular at r = 0. and have to be bound due 
to the u1 term in the KS equation. Far from R they behave asii@ 
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.„.M'^©"'e-'"-''4?(^) . (3) 

where L = \Jh/mujc is the magnetic length and Ll^l is a generahzed Laguerre polynomial^!. 

Once the KS ground state has been obtained, we determine the dipole longitudinal 
response employing linear-response theoryE3. We sketch here how the method is applied. 

For independent electrons in the KS mean field, the variation 5nf^ induced in the spin 
density n„ (a =|, |) by an external field F, whose non-temporal dependence we denote as 
F = J2a fa{^ can be written a^l 

6n^^\r,u;) =Y. j df'x^^l{r,r'-u;)fAr') , (4) 

cr' 

where x^aa' -^S spin density correlation function. In this limit, the frequency u corre- 

sponds to the harmonic time-dependence of the external field F and of the induced 5n^^\ 
Eq. is a 2x2 matrix equation in the two-component Pauli space. In longitudinal re- 
sponse theory, F is diagonal in this space, and its diagonal components are written as a 

vector F = ( \. We consider external dipole (L = 1) fields of the kind 




and FiT) = /(r)e±^M \ ) (5) 




fl 



which cause, respectively, charge and spin density modes of dipole type. To distinguish the 
induced densities in each excitation channel, we label them with an additional superscript 
as 54°'") or 54°'™^- 

The TDLSDFT induced densities are obtained solving the integral equations 

(f, u) = 54°'^) (f, ^) + E / dridf2 xfa, if, n; uj)K„,^, (fi, fa) 5ni^^ (fs, u) , (6) 

(Ti(72 

where either A = n or A = m, and the kernel Kfjfji{r,r') is the electron-hole interaction. 

Equations (^) have been solved as a generalized matrix equation in coordinate 
space. Considering angular decompositions of Xaa' and K^^i of the kind K^^/{r,r') = 
Yl,iKl^\i{r,r')e^^'^^~^'\ it is enough to solve them for I = 1 because only this term couples to 
the external dipole field. 

When a magnetic field is perpendicularly applied to the antidot, the ±1 modes are not 
degenerate and two excitation branches may appear, each having in principle a different 
circular polarization, i.e., carrying an orbital angular momentum ALz = ±1, where Lz is 
that of the gs. In contradistinction with the quantum dot (QD) case, we will see that this 
is not so for antidots, for which the more intense peaks in both branches have the same 
polarization. 

The induced charge or magnetization densities corresponding to density and spin re- 
sponses are given by 6n^^^ = ^raj'^^ + Sn^^^ and 6m^^^ = 6n^^'^^ — ^^j;'^'*. /,From them, the 
dynamical polarizabilities in the density and spin channels are respectively given by 



ann(t^) = j drr f{r)6n^"-\r) 

ammi^)= I drrf{r)5m^^\r) . (7) 
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In this expression Jn'^") has to be understood as the charge density induced by a spin- 
independent operator, and 5m^^^ as the spin density magnetization induced by a spin- 
dependent operator. Within longitudinal response theory, cross-channel induced densities 
such as Sm^""^ or Sn^"^^ may also appeai@. The above expressions hold for each ±1 circular 
polarization of the F operators defined in Eq. (^). Taking into account both possibilities, 
we define a^^\{uj) = aAA{^,^^) + ctAA{—^,^), whose imaginary part is proportional to the 
strength function S^^Ai^) — The peaks appearing in the strength functions 

are dipole charge density (CDE) or spin density excitations (SDE) caused by the external 
field. Analogously, the peaks appearing in the strength function which results from using in 
Eq. the KS induced densities Sn^^'^^ instead of the correlated ones Sn'^J^\ correspond to 
dipole single particle excitations (SPE). 

We have used for the /(r) function entering Eq. (|^) two different choices. One is the 
standard dipole operator /(r) = r, and the other one is /(r) = 1/r. As indicated in Ref. 
the latter choice is inspired in the small-argument expansion of the irregular Bessel function 
Yi(gr). The relevance of this operator in the present context is that it mostly causes edge 
excitations. This will be futher discussed in the next Section. The situation is analogous 
to that faced in the description of surface modes of 3D cavities in metals, see Ref. ^ for a 
thorough discussion. 



III. RESULTS AND DISCUSSION 

As a case of study we have considered an antidot of Us = 0.25 (cto)"^ ^ = '^■^ '^o- 
It roughly corresponds to one of the systems studied by Zhao et ali having Ug = 2.6 x 10^^ 
cm~^ and the same area. We recall that the antidot arrays studied in this reference have a 
very short period, whereas we have taken for R^o defined in Eq. (Q) a mean value of 30 a^, 
larger for large z/'s, and smaller for small i/'s. We have employed B values yielding integer 
filling factors u = 2iTnsC'^ in the range 1 < < 10 which corresponds to 10.8 > B > 1.08 
T. 

Figure shows the particle n(r) and spin magnetization densities m(r) as a function of 
r for = 1 to 6. Away from R, the electron density approaches n^, and either m(r) is zero 
if u is even, or the local spin polarization ^(r) = m{r)/n{r) equals l/u. The number of sp 
orbitals needed to numerically achieve this limit increases with u, making it rather involved 
to obtain high-z/ ground states. For instance, some 2500 occupied sp levels have been used 
to model the u = 10 gs. 

Figure ^ shows the sp energies enia as a function of / for the same configurations as 
in Fig. |I[ In the bulk, enia are arranged into Landau bands characterized by the index 
M = n + — I) /2 and the value of a. The filling factor represents the number of occupied 
Landau bands, each of them labeled as (M, | or J,). It is worth noticing that since we have 
taken B in the positive z direction and the sp orbital angular momentum is written as —I, 
most occupied sp levels have / > 0, and that for a given M the (M, |) band lies below the 
(M, i) one. 

These bands are almost degenerate when u is even because of the smallness of the Zeeman 
energy and because m(r) ~ in the gs. When u is odd, the spin dependence of the exchange- 
correlation energy f„ Vre.po„=ible of the large gap betweek the (M. Ti) bandfl. Dehnig 
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the M-th Landau level as the set of (M, |) and (M, |) bands with M = 0, 1, 2, 
from Fig. |^ that even i> values correspond to completely filled Landau levels configurations, 
and odd u values to configurations in which one Landau level is half-filled. If M is not too 
small, the M-th level is made of sp states having n = M. As in QD's, this rule is violated 
by a few sp states per Landau level at most, see for instance Fig. 1 of Ref. ^ and the u = 6 
panel of Fig. |[ 

The dipole FIR response has been obtained as outlined in Sect. II. Besides evaluating 
the occupied Landau levels, two more empty levels have been calculated, which allows an 
accurate description of modes up to 2 Uc- Figure ^ shows the charge, spin, and single electron 
strength functions corresponding to the dipole operator. All peaks displayed in this figure 
are of cyclotron-like type, i.e., they are excited by the e"*"*^ component of the F operators 
in Eq. ^. This makes possible the transfering of strength among them, an effect that has 
been experimentally observedBH. We remark in passing that the fact that all these peaks 
originate from the same excitation operator makes the sum rule approachBS unsuitable to 
describe their excitation energies. 

Our sign convention implies that at B ^ the gs has a negative total orbital angular 
momentum L^. Thus, these excitations decrease in one unit the absolute value of L^; when- 
ever needed, we shall label cyclotron-like modes with a (— ) sign. Bulk magnetoplasmons 
in QD's also bear this character, whereas edge magnetoplasmons do not (see for instance 
Refs. |26|j3T| ). Indeed, edge magnetoplasmons in QD's are excited by the e~*^ component of 
the dipole field and correspond to excitations which increase in one unit the absolute value 
of L^. Whenever needed, we shall label the anticyclotron-like modes with a (+) sign. The 
cyclotron-like character of edge magnetoplasmons in antidots has its microscopic origin in 
the upwards bending of the Landau levels for small I values, whereas the anticyclotron-like 
character of the edge magnetoplasmons in dots comes from the upwards bending of the 
Landau levels for large I values. 

As in quantum dots, SDK's are weakly redshifted with respect to SPE's because of the at- 
tractive character of the rather small exchange-correlation vertex corrections, whereas CDE's 
are strongly blueshifted because of the repulsive character of the intense direct Coulomb 
interaction. In contradistinction with the QD casell, no clear signature of a collective lon- 
gitudinal spin edge magnetoplasmon appears at small B values: the edge spin mode is 
rather fragmented and weak. We recall that at full magnetizacion longitudinal spin and 
density responses coincideil and only spin modes of transverse type, not studied here, can 
be excitedc3. 

Another interesting feature of longitudinal response theory is that spin and charge density 
modes are generally coupled if the gs has a large spin magnetization, as it happens for 
instance at = 3. This means that spin dependent probes may excite charge density modes. 
Conversely, spin independent probes may excite spin density modes. This coupling has been 
experimentally detected in QD's using Raman spectroscopy!!, and has been theoretically 
addressed^ within TDLSDFT. For an antidot the effect is more marked, since ^ reaches l/u 
instead of a smaller effective value as it happens in QD. This can be seen in Fig. ^. 

It is worthwhile to compare antidot with dot results. For this sake, we show in Fig. 
^ the dipole charge density response of a QD made of = 210 electrons confined by the 
potential created by a jellium disk of radius R ~ 16.7 Oq. Outside the edge region, this 
system has a surface density of ~ 0.24(aQ)~^, similar to our case-of study antidot, and it 
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has been thoroughly studied in Ref. For this dot the 1 < z/ < 10 range corresponds to 
10.3 >B> 1.03 T. 

The energies of the more intense CDE's of the antidot and of the = 210 QD have 
been drawn as a function of B in Figs. |^ and respectively (a more complete version 
of Fig. 1^ can be found in Ref. |2^). Also drawn are the cyclotron frequency Uc, and 2uJc 
and 3uJc as well. These figures show that the interaction of the bulk magnetoplasmon 
with the harmonics of the cyclotron resonance ntOc, resembling the Berstein modesEa, which 
causes the splitting of the magnetoplasmon near 2u)c and 3uJc, is well reproduced by the 
calculations. Alongside the groups of peaks corresponding to CDE's we have indicated the 
change in radial quantum number n of the sp wave functions involved in the transition. This 
change An is unambiguosly identified in TDLSDFT calculations introducing energy cutoffs 
in the electron-hole energy denominators Ae entering the definition of the free correlation 
functionii x^aa'- limiting the electron-hole pairs to those having Ae < Uc, only An = peaks 
apear in the response, whereas limiting them to those with Ae < 2uJc, no An = 2 peaks 
appear. These cyclotron-like An = 2 peaks have been experimentally detectedi. 

The correlations introduced by TDLSDFT have a dramatic effect on the FIR response 
of the system which might not always be properly recognized. Not only they appreciably 
shift the CDE energies with respect to the SPE's but also affect their intensity, causing that 
peaks having a sizeable strength in the SPE response are nearly washed out, whereas others 
hardly perceived in the SPE response adquire a large strength. To illustrate it we show in 
Fig. § the single particle and charge density strength functions for z/ = 7 in a logarithmic 
scale. It can be noticed that in the SPE strength the 2 Uc peak, which has an intensity six 
orders of magnitude smaller than the Uc peak and cannot be seen at the scale of Fig. |^, 
emerges as a very collective peak in the CDE strength, of intensity similar to that of the 
bulk An = 1 peak. Similarly, the weak edge SPE's loose their dimy intensity, transferring 
it to the very collective An = edge CDE. Of course, in actual TDLSDFT calculations all 
these modes are reciprocally influenced, but the gross effect is as described. 

TDLSDFT does reproduce the oscillation of the edge magnetoplasmon energy as a func- 
tion of filling factoiH. Figure ^ shows the edge magnetoplasmon energy as a function of u 
for the antidot and the = 210 dot. The oscillations and the out-of-phase effectlll are tiny 
effects better seen plotting the second differences of the edge magnetoplasmon energy: 

A2cu(z/) =tu(z/ + l)-2cu(z/)+cj(z/-l) . (8) 

The energy oscillations are clearly seen in the Fig. ^ insert, as well as the experimental 
finding that maxima in the edge magnetoplasmon energies correspond to half-filled Landau 
levels in antidots, and to filled Landau levels in dots. 

^From Figs. ^ and |^ we conclude that for an isolate antidot, TDLSDFT is unable to 
produce the anticrossing of the edge and the cyclotron modes at small B values which is 
observed in short-period antidot arrays. Yet, the behavior of the antidot and dot An = 
edge branch near the crossing is quite distinct, manifesting a different curvature. The low- 
B behavior of the edge brach experimentally found in short-period antidot arrays is likely 
influenced by neighboring antidots. Notice that C ~ 1.9 Og at z^ = 6, which is roughly half 
the effective surface width of the antidot lattice period of Ref. 0. 

Figures ^ and ^ also show that the TDLSDFT antidot response is very fragmented at 
low B. This renders tricky to distinguish bulk from edge charge modes. Notice that at 
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z/ = 10 the high energy peak in Fig. ^ is a An = 2 mode, whereas the low energy peak is a 
mixture of An = and An = 1 peaks. The situation worsens increasing u, which is beyond 
our capabihties. For instance, to have results at S ~ 0.5 T one would need to carry out 
calculations for M ~ 20. 

A possible way to partially disentangle these modes is to study the charge density re- 
sponse to the /(r) = 1/r operator. Its r dependence makes this very ineffective in exciting 
cyclotron-like modes while it is well suited to describe edge oscillations. This can be seeing 
in Fig. where we have plotted the CDE's caused by r and 1/r for z/ = 10 to 7. In 
the bottom panel, the intensity of the 1/r peaks has been multiply by a factor of 50,000 to 
allow for a sensible comparison with the r peaks. It can then be seen in that panel that 
both operators excite the same AL^ = +1 peaks, but with much different intensity, and 
that indeed, the low energy peak clearly visible in the response is an edge peak. 

Finally, we have plotted in the top panel of that figure the CDE caused by the e~*^ 
component of the F operators in Eq. ^. These are anticyclotron-like excitations of very 
weak intensity: to have them at the same scale as in the bottom panel, one should divide 
the scale of the r modes by a factor of 40, and the scale of the 1/r modes by a factor of 
4, i.e., the 1/r operator excites the cyclotron-like and anticyclotron-like main peaks with 
roughly equal intensity. In agreement with the experiment (see especially Fig. 1 of Ref. 
it can be seen on the one hand that their intensity quickly decrease as B increases, and on 
the other hand that their energies are very similar to these of the cyclotron-like modes. This 
makes these modes experimentaly accessible only using circular polarized radiation. 
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FIGURES 





FIG. 1. Electron (solid lines) and spin magnetization (dashed lines) densities as a function of 
r for an antidot of R = 7.5 and = 0.25 (ag)"^ and several filling factors. Also shown is the 
jellium density (dotted lines). 
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FIG. 2. Single particle energies enia as a function of / for the configurations of Fig. ||. The 
horizontal line represents the electron chemical potential. Full upright triangles represent (M, |) 
bands, and empty downright triangles represent (M, |) bands. 
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FIG. 3. Charge, spin, and single particle antidot dipole strength functions in arbitrary units. 
All displayed modes are cyclotron-like polarized (— ). Notice the different energy scale of the u = 1 
responses. 
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FIG. 4. Charge and spin antidot dipole strength functions at v = 3 in arbitrary units illus- 
trating the coupling between the strengths at high spin magnetizations. 
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FIG. 5. Charge and single particle dipole stregth functions in arbitrary units corresponding to 
the A'^ = 210 dot described in the text. Cyclotron-like polarized (— ) bulk modes and anticyclotron 
polarized (+) edge modes appear in the responses. Notice the different energy scale of the v = 1 
responses. 
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FIG. 6. Energies of the more intense antidot CDE's (solid circles) and SDE's (open circles) as 
a function of B. Prom left to right, the calculated points correspond to = 10 to 2. 
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FIG. 8. Antidot CDE (solid line) and SPE (dashed line) strength functions at = 7 in a 
logarithmic scale. 
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FIG. 9. Edge magnetoplasmon energies of the antidot (solid circles) and of the N = 210 dot 
(open circles) as a function of filling factor. The insert displays the second differences A2a;. The 
lines have been drawn to guide the eye. 
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FIG. 10. Antidot charge density strength functions in arbitrary units for = 10 to 7 caused 
by r (solid hues) and 1/r (dashed lines) fields. Top panel, anticyclotron-like polarized (+) modes. 
Bottom panel, cyclotron-Hke polarized (— ) modes. 



20 



